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Abstract. Quick chemical equilibration times of hadrons within a hadron gas are 
explained dynamically using Hagedorn states, which drive particles into equilibrium 
close to the critical temperature. Within this scheme master equations are employed 
for the chemical equilibration of various hadronic particles like (strange) baryon and 
antibaryons. A comparison of the Hagedorn model to recent lattice results is made 
and it is found that for both Tc =176 MeV and Tc=196 MeV, the hadrons can 
reach chemical equilibrium almost immediately, well before the chemical freeze-out 
temperatures found in thermal fits for a hadron gas without Hagedorn states. 



1. Introduction and Motivation 

As two heavy ions collide color neutral clusters are formed within which the number of 
particles per cluster increase. The clusters become so dense and begin to overlap such 
that it impossible to distinguish quarks from one cluster from that in another i.e. a 
percolation transition. The critical density for this is about e ~ 1 GeV / fm? . Following 
the phase transition into Quark Gluon Plasma the interactions are dominated by quarks 
and gluons. Through gluon fusions, strange quarks can easily be reproduced. Eventually 
the QGP cools back into hadrons where the particle yields and ratios are measured. 

If one only considers binary collisions, which react too slowly for strange particles 
to reach chemical equilibrium within the hadron gas phase, then it is clear that strange 
particle yields can only be explained through gluon fusion within QGP P and that 
the hadrons must exist QGP already in full chemical equilibrium [2|. However, multi- 
mesonic collisions mx ^ XX have been demonstrated to reach chemical equilibration 
for varoius (strange) antibaryons quickly at SPS [3] , although they are still not enough 
to explain the particle yields of exotic antibaryons at the higher energies at RHIC 
[U [5]. In order to circumvent such longer time scales ~ 10 fm/c for a situation of 
a nearly baryon-free system with nearly as much antibaryons as baryons, it was then 
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suggested by Braun-Munzinger, Stachel and Wetterich [6] that near Tc there exists an 
extra large particle density overpopulated with pions and kaons, which then drive the 
baryons/anti-baryons into equilibrium by exactly such multi-mesonic collisions. But it 
is not clear how and why this overpopulation of pions and kaons should appear, and how 
the subsequent population of (anti-)baryons would follow in accordance with a standard 
statistical hadron model: According to the mass action law the overpopulated matter 
of pions will result in an overpopulation of (anti-)baryons. For such a large number of 
(anti-)baryons it is difficult to get rid of them quickly enough in order to reach standard 
hadron equilibrium values before the chemical freeze-out [7]. 

Rather, understanding the rapid chemical equilibration is possible using Hagedorn 
states, heavy resonances that drive similar and more multi-hadronic reactions close to 
Tc, as shown in [71 |Hl El [IH] ■ Close to the matter is then a strongly interacting mixture 
of standard hadrons and such resonances. Using the Hagedorn states as potential and 
highly unstable catalysts, the standard hadrons can be populated reactions: 

nn^ HS ^ n'n + XX (1) 

where XX can be substituted with pp , KK , AA, or QQ. The large masses of the 
decaying Hagedorn states open up the phase space for multi-particle decays. 

In this note we will compare the particle ratios obtained by using reactions driven 
by Hagedorn states and those of the experiments at RHIC. We find that both strange 
and non-strange particles match the experimental data well within the error bars. 
Furthermore, we are able to make estimates for the chemical equilibration time and 
find that they are very short, which implies that the hadrons can easily reach chemical 
equilibrium within an expanding, hadronic fireball and that hadrons do not need to be 
"born" into chemical equilibrium [HI HJ [10] . Hagedorn states thus provide a microscopic 
basis for understanding hadronisation of deconfined matter to all hadronic particles. 

Before starting with the details, we emphasize that Hagedorn states have become 
quite popular to understand the physics of strongly interacting matter close to the 
critical temperature: Hagedorn states have been shown to contribute to the physical 
description of a hadron gas close to T^. The inclusion of Hagedorn states leads to a low 
7]/s in the hadron gas phase [H], which nears the string theory bound rj/s = 1/{Ait) 
. Calculations of the trace anomaly including Hagedorn states also fits recent lattice 
results well and correctly describe the minimum of the speed of sound squared, c^, near 
the phase transition found on the lattice [TT] . Estimates for the bulk viscosity including 
Hagedorn states in the hadron gas phase indicate that the bulk viscosity, (/s, increases 
near Tc [Hj. We also remark here that Hagedorn states can also explain the phase 
transition above the critical temperature and, depending on the intrinsic parameters, 
the order of the phase transition [12]. Finally, it has been shown [13] that Hagedorn 
states provide a better fit within a thermal model to the hadron yield particle ratios. 
Additionally, Hagedorn states provide a mechanism to relate and Tchem, which then 
leads to the suggestion that a lower critical temperature could possibly be preferred, 
according to the thermal fits [TB] . 
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2. Setup 

The basis of the Hagedorn spectrum is that there is an exponential mass increase along 
with a prefactor i.e. the mass spectrum has the form: f{m) !^ exp™/^^ [H]. The 
exponential mass spectrum drives open the phase space, which allows for multi-mesonic 
decays to dominate close to Tc (we assume Th ~ Tc). We use the form 

r^^ A ^ 

p= ^e^H dm. (2) 

J Mo [nfi + m^] 4 

where Mq = 2 GeV and = 0.5 GeV. We consider two different different lattice 
results for T^.: T^. = 196 MeV dn] (the corresponding fit to the trace anomaly is then 
A = 0.5GeV^/^, M = 12 GeV, and B = (340MeV^)^), which uses an almost physical 
pion mass, and = 176 MeV [T7j (the corresponding fit to the energy density leads to 
A = O.lGel/3/^ M = 12 GeV, and B = (300Me\/)^). Both are shown and discussed 
[To] . Furthermore, we need to take into account the repulsive interactions and, thus, use 
volume corrections [101 dH HH], which ensure that the our model is thermodynamically 
consistent. Note that i? is a free parameter based upon the idea of the MIT bag constant. 

We need to consider the back reactions of multiple particles combining to form 
a Hagedorn state in order to preserve detailed balance. The rate equations for the 
Hagedorn resonances iVj, pions iV^, and the XX pair Nxx, respectively, are given by 
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The decay widths for the i*^ resonance are Fj ,^ and T^xxi the branching ratio is Bi n 
(see below), and the average number of pions that each resonance will decay into is (nj). 
The equilibrium values N^'^ are both temperature and chemical potential dependent. 
However, here we set ph = 0. Additionally, a discrete spectrum of Hagedorn states is 
considered, which is separated into mass bins of 100 MeV. 

The branching ratios, Bi^n, are the probability that the i^^ Hagedorn state will 
decay into n pions where ^„ Bi^n = 1 must always hold. We assume the branching 
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ratios follow a Gaussian distribution for the reaction HS o rnr 

1 



5,,^— =e , (4) 



which has its peak centered at (n^) ^3 — 34 and the width of the distribution is 
^ 0.8 — 510 (see (ID])- For the average number of pions when a XX pair is present, 
we again refer to the micro-canonical model in [71 [19] and find 

(ui,,) = ( ^ ) (0.3 + OAmi) ^2-7. (5) 



where rrii is in GeV. In this paper we do not consider a distribution but rather 
only the average number of pions when a XX pair is present. We assume that 
= {'^i,p) = {ni,k) = {^i,A) = ioT wheu a kaon anti-kaon pair, AA, or Qfl 

pair is present. 

The decays widths are defined as follows (see [TU]): 

Ti = O.lSmi - 0.0584 = 250 MeV to 1800 MeV 
Ti^pp = 3 MeV to 1000 MeV 
T^ j^j^ = 50 MeV to 1700 MeV 
r^^^^ = 3 MeV to 250 MeV 
Ti^^^ = .01 MeV to 4 MeV 

Ti,^ =Ti-Ti^xx- (6) 

Fj is a linear fit extrapolated from the data in . It is then seperated into two parts, 
one for the reaction HS -H- nn i.e. Fj^^r and one for the reaction HS ^ nvr + XX i.e. 
T^ xx- The decay width T^^xx is found my multiplying {X), which is the average X that 
a Hagedorn state of mass m will decay into, that is found from both microcanonical 
[ini [7] and canonical models |T0]. The large masses open up the phase space for such 
more special multi-particle decays. A detailed explanation is found in [TU] . 

The equihbrium values are found using a statistical model [21], which includes 104 
light or strange particles from the the PDG [20] • As in [21], we also consider the effects 
of feeding for pions, protons, kaons, and lambdas. Additionally, throughout this paper 
our initial conditions are the various fugacities at t^ (at the point of the phase transition 
into the hadron gas phase) a = A^(to) , A = Kito) ,and0 = Xxxi^o) which are chosen 
by holding the contribution to the total entropy from the Hagedorn states and pions 
constant i.e. SHad{To,a)V{to) + shs{To, l3i)V{to) = SHad+HsiTo)V{to) = const and the 
corresponding initial condition configurations we choose later can be seen in Tab. |l| (for 
further discussion see [TO]). 



3. Results: Expansion 



In order to include the cooling of the fireball we need to find a relationship between the 
temperature and the time, i.e., T{t). To do this we apply a Bjorken expansion for which 
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the total entropy is held constant 

const. = s{T)V{t) ~ ^ I ^dy. (7) 

where s(T) is the entropy density of the hadron gas with volume corrections. The total 
number of pions in the 5% most central collisions, can be found from experimental 

results in [22]. Thus, our total pion number is J^i^tt' = I~o5 ^^^V = While for 



dy 

a gas of non- interacting Bose gas of massless pions St^/N^ = 3.6, we do have a mass 
for a our pions, so we must adjust S^^/Nt^ accordingly. In (2^ it was shown that when 
the pions have a mass the ratio changes and, therefore, the entropy per pion is close to 
Stj/Nj, ^ 5.5, which is what we use here. 

The effective volume at mid-rapidity can be parametrized as a function of time. 
We do this by using a Bjorken expansion and including accelerating radial flow. The 
volume term is then 

V{t) = 71 at (^ro + vo{t - to) + ^ao{t - tof^ (8) 

where the initial radius is ro(to) = 7.1 fm for Th = 196 and the corresponding 
t'^^^^ ^ 2/m/c. For Th = 176 we allow for a longer expansion before the hadron gas 
phase is reached and, thus, calculate the appropriate t^o'^^^ from the expansion starting 
at Th = 196, which is t^^'^^^ ^ 4/m/c. We use Vo = 0.5 and Qq = 0.025 (see |10j). 

Because the volume expansion depends on the entropy and the Hagedorn resonances 
contribute strongly to the entropy only close to the critical temperature (see [ID]), the 
effects of the Hagedorn states must be taken into account with calculating the total 
particle yields otherwise the yields do not increase with the temperature (see [10] for 
further discussion). Therefore, one has to include the potential contribution of the 
Hagedorn resonances to the pions as in the case of standard hadronic resonances, e.g. a 
p-meson decays dominantly into two pions and, thus, accounts for them by a factor two. 
Including Hagedorn states, we arrive at our effective number of pions and XX pairs 

Kxx = N^ + J2Ni [(1 - (X,)) (n,) + (X,)(n,.)] 

i 

Nxx =iV^x + E^^(^^) (9) 

i 

because Hagedorn states also contribute strongly to the XX pairs close to Tc. 

Along with the expansion we also must solve these rate equations, Eq. ([3]), 
numerically . We start with various initial conditions as seen in table II and the initial 
temperature is the respective critical temperature and we stop at T = 110 MeV. The 
black solid line in each graph is the chemical equilibrium abundances and the colored 
lines are the various initial conditions listed in Tab. [T} In order to save space, the results 
are only the expansion of the A's and fi's results are only shown for Th = 196 MeV. 
However, Th = 176 MeV and further results are shown and discussed in [lOj. However, 
the end particle ratios are all shown in Fig. |2} Note that in all the following figures the 
effective numbers are shown so that the contribution of the Hagedorn states is included. 



Thermalization through Hagedorn states - the importance of multiparticle collisions 6 







ft. — \ .(fA 


V — ^xxy^o) 




1 


1 





IC2 


1 


1 


0.5 


ic^ 


1.1 


0.5 





Id 


0.95 


1.2 






Table 1. Initial condition configurations. 
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Figure 1. Results for A's and fJ's with various initial conditions for Th = 196 MeV. 



We can also observe the affects of the expansion on the AA pairs as seen in Fig. 
[1] for the reaction nn o HS ^ nn + AA. They reach the experimental values almost 
immediately. One can see that the chemical equilibration time does depend slightly 
on our choice of Pi, i.e., a larger /3j means a quicker chemical equilibration time. The 
one exception is for an underpopulation < 1) of Hagedorn states, which reaches 
chemical equilibrium by T = 170 for = 196 MeV. Additionally, when the AA pairs 
start at about half their chemical equilibrium values, it only helps the AA pairs to reach 
equilibrium at a slightly higher temperature (on the order of a couple of MeV). 

We used our model to investigate the possibility of fi's produced through [7]: 

HS ^m + x 

HS (sssqqq) ^n + B + X 

HSb{sss) ^n + X. (10) 

the first we implemented into our model to obtain Fig. [T] where we are impressively able 
to populate the fifi's. However, this scenario is slightly more delicate because F, by 
50%, the total production of Q is not sufficient up to 25 %, to meet the experimental 
yield (the other ratios are not significantly affected by the change in Fj). 

A summary graph of all our results is shown in Fig. [2j The gray error bars cover the 
range of error for the experimental data points from STAR and PHENIX. The points 
show the range in values for the initial conditions at a final expansion point with a 
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Figure 2. Plot of the various ratios including all initial conditions defined in Tab. [T] 
The points show the ratios at T = 110 MeV for the various initial conditions (circles 
are for Th = 176 MeV and diamonds are for Th — 196 MeV). The experimental results 
for STAR and PHENIX are shown by the brown error bars. 



temperature T = 110 MeV. We see in our graph that our freezeout results match the 
experiments well and the initial conditions have little effect on the ratios, which implies 
that information from the QGP regarding multiplicities is washed out due to the rapid 
dynamics of Hagedorn states. A smaller /3j slows the equilibrium time slightly. However, 
as seen in Fig.[2]it still fits within the experimental values. Furthermore, in [8] we showed 
the the initial condition play almost no roll whatsoever in K/tt^ and {B + B)/tt~^. Thus, 
strengthening our argument that the dynamics are washed out following the QGP. 



4. Conclusion 



To conclude we have found that the Hagedorn states provide a mechanism for quick 
chemical equilibration times. Our model gives chemical equilibration times on the 
order of At ^ 1 — 3^. Furthermore, the particle ratios obtained from decays of 
Hagedorn states match the experimental values at RHIC very well, which leads to the 
conclusion that hadrons do not need to be born in chemical equilibrium. Rather a 
scenario of hadrons that reach chemical freeze-out shortly after the critical temperature 
due to multi-mesonic reactions driven by Hagedorn states, is entirely plausible. We 
have shown that both strange (A's and K's) and non-strange (vr's and p's) hadrons 
can reach chemical equilibration by T = 160 MeV. Thus, it would be interesting to 
implement Hagedorn states into a transport approach such as UrQMD |2l] Moreoever, 
even multi-strange baryons such as fi's can reach chemical equilibrium in such a scenario. 
However, unlike the other hadrons, the final abundancy of Qs is slightly dependent on 
the decay width of the Hagedorn states (they cannot reach full chemical equilibrium if 
the decay width is 50% of its current value) and, thus, should be further investigated 
using strange and/or baryonic Hagedorn states J7]. Although they are quite exotic, it 
is still possible that they might occur. Still, our work indicates that the population 
and repopulation of potential Hagedorn states close to phase boundary can be the key 
source for a dynamical understanding of generating and chemically equilibrating the 
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standard and measured hadrons. Hagedorn states thus can provide a microscopic basis 
for understanding hadronisation of deconfined matter. 
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